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INTRODUCTION 

Studies of the radiation balance of the Earth-atmosphere system and its components provide a 
means for investigating the energy exchange between the planet Earth and space and the radiative 
energy available for driving the atmospheric circulation. They further provide a check of the energy 
balance of atmospheric circulation models at their upper boundary and make possible a variety of 
climatological investigations and radiation budget studies at various locations of the Earth, which help 
quite considerably to understand the nature of the atmosphere’s general circulation. 

A fundamental requirement for such studies is the existence of complete global observations at a 
proper spatial resolution having high accuracy and extending over a period of several years. Previous 
satellite observations made with relatively high resolution instruments, such as the medium resolution 
infrared radiometer (MRIR) flown on various TIROS satellites and on Nimbus 2, did not meet, for 
several reasons, this requirement in all respects. The measurements from Nimbus 3 extended over a 
continuous period of nearly 5 months plus two more semimonthly periods permitting estimates of the 
annual global budget from scanning radiometer measurements for the first time. 

It is the purpose of this report to present in detail a description of the evaluation methods used, 
a discussion of the basic assumptions inherent in them, and various illustrations of possible sources of 
systematic errors. Results are presented only to illustrate the annual cycle of the measured radiation 
budget. Broader aspects of these studies will be discussed at greater length in subsequent papers. A few 
preliminary results from the Nimbus 3 experiment have already been published (Raschke et al., 1971 ; 
Vonder Haar and Raschke, 1972; and Vonder Haar et al., 1972. 



DEFINITIONS 


In all discussions of this report the Earth-atmosphere system at each location X, 0 (geographic 
longitude and latitude) is considered to be an infinitesimally thin layer on top of a spherical Earth with 
mean radius R Q = 6371 km. Daily averages of the net flux of radiative energy across a horizontal area 
element at the top of this layer, in the spectral range between 0.2 and 200 jum, can be computed as 
follows: 

Q{\, 0, t) = H s (k, 0, t ) - ,4(X, 0, t) ■ H s (K, 0, t) - W L (X, 0, t) (1) 

The quantity Q will be called hereafter the radiation balance of the Earth-atmosphere system. Its com- 
ponents are— 

i/ s (X, 0, t) Daily average of the irradiance (for its definition, see Nicodemus, 1963) of the 
extraterrestrial, electromagnetic solar radiation arriving at the Earth, values of 
which have been computed here assuming the solar constant to be S Q = 1.95 cal 
cm -2 ' min - 11 (Drummond, 1970). The flux of particle radiation is some orders of 
magnitude smaller. Units are calories per square centimeter per minute. 

A(X, 0, t) Daily average of the albedo of the Earth-atmosphere system in the spectral range 
between 0.2 and 3.0 gra. Units are dimensionless. 

W L (X, 0, t ) Daily average of the outgoing longwave radiation or of the radiant exitance of 

the Earth-atmosphere system to space in the spectral region between about 4.0 
and 200 jam. Units are calories per square centimeter per minute. 

t Time, referring to the day of the year. 

The second term in equation (1) is the flux density of shortwave radiation reflected and scattered 
back to space by the Earth’s surface, clouds, air molecules, and dust within the atmosphere. Thus, the 
difference between the first and second term of equation (1) describes the amount of radiative energy 
absorbed per time unit in a column of the Earth-atmosphere system. The outgoing longwave radiation 
W L originates from the Earth’s surface and clouds, particulates, and the optically active gases within 
the atmosphere. Both reflection and emission mostly arise from processes in the troposphere and strato- 
sphere, while only very small fractions of the incident radiation are absorbed in the upper atmospheric 
layers above a level of about 60 km. Also the emission from those upper layers to space is only a small 
fraction of W L . Thus the radiation balance Q, as determined from equation (1), describes primarily the 
exchange of electromagnetic radiation between space and the Earth-atmosphere system below a level 
of about 50 to 60 km above the Earth’s surface. 

A relative comparison of the contributions of various components of the Earth-atmosphere sys- 
tem to the radiation balance can be shown by reference to London’s (1957) calculations. Table 1 
shows that more than 80 percent of reflected solar and emitted thermal radiation originates from 
clouds and tropospheric air. 

All quantities in equation ( 1 ) except the albedo are radiant flux densities at a horizontal element 
of area. The general notation is 
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Table 1.— Annual Radiation 

Budget of 

the Earth- 

Atmosphere System 


[after London, 1957] 



cal cm -2 
min -1 

percent 

Insolation at top of the atmosphere 

0.5 

100 

Reflection and scattering to space: 

0.176 

100 

From atmosphere 

0.034 

19 

From clouds 

0.121 

69 

From Earth’s surface 

0.021 

12 

Longwave radiation lost to space: 

0.324 

100 

From troposphere 

0.283 

88 

From stratosphere 

0.014 

4 

From Earth’s surface 

0.027 

8 


n rr/2 

N(6, \p) cos 9 sin 9 dd d\p (2) 

. 

where N(6, i//) is the radiance impinging on (or emerging from) that surface element from directions 
6 and \p. 6 and \ {/ are the zenith and azimuth angles, respectively, from the principal plane, which is 
the plane including the Sun, the Earth’s center, and the location of observation. (See app. A; fig. A-l.) 

The albedo of the Earth-atmosphere system is defined in this report as the ratio between daily 
averages of outgoing (reflected and scattered) and incoming solar radiation in the spectral range between 
0.2 and 3.0 jam. This quantity corresponds physically to the directional reflectance of a surface or the 
Earth-atmosphere system, which is the ratio between outgoing and incoming radiant energy per area 
and time unit, but at the time of observation. The adjective “directional” refers here to the zenith 
angle of an illuminating radiation source, such as the Sun. If only the radiance N within a narrow field 
of view is observed, then one can determine the reflectivity of a surface into the observational direction 
only, the bidirectional reflectance 


p(9, \p, ?) = 


H s cos f 


sr 


.- 1 


(3) 


where f = zenith angle of the source of illumination, and 9 and i// are the zenith and azimuth angles of 
observation. H s is required here to be a plane wave or parallel radiation, which can be considered to be 
the case for the Sun’s radiation above the atmosphere. From equation (3) the directional reflectance r 
of a surface of the Earth-atmosphere system at a specific location can be obtained by integration ac- 
cording to equation (2): 


*£) = 



p(9, i p, f) cos 9 sin 9 d9 d\p 


(4) 
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Over a complete diffusely reflecting surface, p is independent of 6 and \p, thus its directional reflect- 
ance is 

r d (£) = irp($) (5) 

This quantity will be used in the section of this report entitled “Data Accuracy and Precision” 
for preliminary checks of measurements of reflected solar radiation from Nimbus 3. There it will be 
called simply reflectance; it is not the same as the albedo A of equation (1 ). 

AVAILABLE DATA 

Nimbus 3 Medium-Resolution Infrared Radiometer Experiment 

The satellite Nimbus 3 was launched on April 14, 1969, into a retrograde, Sun-synchronous, and 
nearly circular orbit. Its orbital characteristics were as follows: 

Perigee: 1097 km 

Apogee: 1143 km 

Inclination: 99.9° 

Orbital period: 107.3 min 

Nodal time: 11:30 a.m. (local time) 

These characteristics specified longitudinal displacements of 26.8° between each equatorial node. 

Thus, during a 24-hr period, measurements could be obtained over the entire globe, once under day- 
light and once under nighttime conditions, at about local noon or midnight, in tropical and subtropi- 
cal regions (Nimbus Project Staff, 1969a, b). In addition to other experiments for meteorological and 
physical purposes, a cross-track scanning five-channel MRIR of the type flown on Nimbus 2 and sev- 
eral TIROS satellites (McCulloch, 1969) was flown on this satellite. 

Data were available for the entire period between April 16 arid August 15, 1969. Malfunction of 
one of the tape recorders allowed recording of only two further semimonthly periods, which were se- 
lected specifically to complete an entire annual cycle (October 3 to 17, 1969, and January 21 to Feb- 
ruary 3, 1970). This set of available data covers one entire season and some adjacent weeks. The two 
further periods of data, of course, do not cover the entire seasons from which they were taken, but 
they allow at least an estimate of the annual cycle and of annual averages of the radiation balance and 
its components. Data taken between May 16 and July 30, 1969, allow detailed comparative studies to 
be made with previous measurements by Nimbus 2 in 1966 (Raschke and Bandeen, 1970). 

The five channels of the MRIR experiment were sensitive to radiation in five different spectral 
intervals, whose total and half-power wavelengths are listed in table 2 together with some descriptive 
characteristics of the radiation received. Curves of the spectral sensitivity of each channel are shown 
in figures 1 and 2. 

The spectral response of channel 5 covers almost the entire spectral range of incident solar radia- 
tion (fig. 1); however, in the visible and near infrared, channel 5 is only half as sensitive as it is at 
wavelengths greater than 1.5 pm. Thus, spectral reflection properties of observed areas in the infrared 
will be slightly overweighted in the average reflectivity value. This may affect especially the reflectance 
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Table 2.— Spectral Response of the 5-Channel MRIR Flown in Nimbus 3 


Channel 

Total wavelength interval, (im 
(Half-power wavelength interval, /am) 

Notes 

1 

6.0 to 7.0 (6.35 to 6.72) 

Radiation from upper 
tropospheric water vapor 
and cirrus clouds 

. 2 

9.1 to 12.1 (10.1 to 11.2) 

Radiation related to surface 
temperature 

3 

14.0 to 16.3(14.5 to 15.8) 

Radiation related to mean 
temperature of lower 
stratosphere 

4 

20.2 to 23.9 (20.8 to 23.2) 

Radiation from lower 
tropospheric water vapor 

5 

0.2 to 4.8 (0.45 to 3.9) 

Radiation from reflectance 
of surfaces and scattered 
radiation 



0.0 .5 1.0 1 .5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

WAVELENGTH (ji m) 

Figure 1.— Spectral sensitivity of channel 5 and the extraterrestrial spectral irradiance of the Sun. 


computed from measurements taken over such areas, where the spectral reflectivity at wavelengths 
shorter than about 1.0 gm is completely (higher or lower) different from that at longer wavelengths. 
Snow surfaces, whose albedo is very high in the visible b>ut low in the near infrared may, then, be 
slightly underestimated in their reflection properties. Accurate corrections for such errors might be 
possible only if spectral measurements of reflected solar radiation taken at very high altitude and 
different solar zenith angles and over various areas were available. 
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‘ 5.0 7.5 10.0 1 2.5 15.0 1 7.5 20.0 22.5 25.0 

WAVELENGTH (|>n) 

Figure 2.-Spectral sensitivity of channels 1 to 4 and the spectral radiant exitance of a blackbody 

at temperatures of 200, 250, and 300 K. 


Table 3.— Average and Nimbus 3 “Filtered” Reflectances of Artificial 
and Natural Surfaces for Solar Radiation of the Spectral Irradiance 


Surface 

Reflectance 

0.2 to 0.8 pm 

0.8 to 4.8 pm 

Average 

Filtered 

1 

1.0 

0.0 


0.46 

2 

1.0 

0.5 


.74 

3 

.5 

1.0 

.72 

.78 

4 

.0 

1.0 

.44 

.54 

Snow and atmosphere near Cleveland 

.51 

.48 

Snow and atmosphere near Bear Lake, Utah 

.49 



.47 


An estimate of possible errors due to the spectral sensitivity of channel 5 in the whole spectral 
range between 0.2 and 4.8 pm is summarized in table 3. There, average diffuse reflectances of various 
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Figure 3.— Spectral reflectance of snow and cloudfree atmosphere obtained from airborne measurements at different 

altitudes above snow-covered Bear Lake, Utah (Hovis, 1971). 


artificial (cases 1 to 4) and snow surfaces are compared with those that would be obtained from chan- 
nel 5 measurements. In both calculations the spectral diffuse reflectances ( r d , eq. (5)) were weighted 
with the spectral solar irradiance derived recently by Thekaekara (1970). The discussed measurements 
over snow are values of the reflected solar radiation obtained by Hovis et al. (1971) from an altitude 
of 9 to 12 km. Spectral reflectances derived from them are shown in figure 3. They are not corrected 
for the anisotropic nature of the reflection properties of the Earth-atmosphere system. 

The artificial surfaces 1 and 2 were chosen to have a drastically higher reflectance for radiation of 
short wavelengths. The filtered measurements will always slightly underestimate the mean reflectance 
in these cases, while in reverse cases (higher reflectance at longer wavelengths, surfaces 3 and 4), they 
will always overestimate the mean reflectance. 

The spectrum of reflected solar radiation, as it could be observed at the top of the atmosphere or 
from a very-high-flying airplane, shows over almost all types of surfaces higher radiances in the blue 
than in the infrared regions because of the wavelength dependences of Rayleigh and dust scattering in 
the atmosphere. Thus, atmospheric interference will smooth out such drastic contrasts in the spectral 
reflectance of surfaces. Further, the spectral irradiance of the Sun as shown in figure 1 is weighted 
clearly in favor of shortwave components of the spectrum. One should, therefore, obtain from all 
channel 5 measurements a near representative, or, in some cases only very slightly underestimated, 
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Figure 4.— Trend of reflectance over the Sahara Desert. 


average reflectance of the Earth-atmosphere sytem. Again, exact checks of this assumption are possi- 
ble only with experimental data from very-high-flying airplanes or from special satellite experiments. 

None of the four channels of MRIR sensitive in the infrared covered the entire spectral range 
from about 4 to 200 pm, where almost all infrared radiative energy is emitted to space. Each of them 
covered only a very narrow spectral range located at different areas of the spectrum, as shown in fig- 
ure 2. All four simultaneous radiance measurements were combined to determine the total outgoing 
radiance from regression equations that were obtained from model calculations, as will be described 
in detail. 

Data Accuracy and Precision 

The random error of radiance measurements of all five MRIR channels is placed by the experi- 
menters to be not larger than 1 to 2 percent. This number also includes errors caused by other proce- 
dures, such as the digitization of the data. 

The measurements of all four infrared channels were calibrated in space after each scan by com- 
parison with an internal calibration source (or an equivalent pulse) and against radiation from space. 

No relative calibrations were possible for measurements of the reflected solar radiance (channel 5). 
Thus, it was decided to check their trend over a bright and almost cloudfree area whose reflection 
properties were assumed to be constant over the entire year. Such areas are the deserts of North 
Africa and Arabia. A comparison of single measurements as well as averages within grid fields of about 
500 by 250 km 2 (fig. 4) showed no systematic decrease or increase of reflectances during the entire 
period of available observations. These areas were identified to be cloudfree or nearly cloudfree accord- 
ing to simultaneous channel 2 observations. 

No checks of changes in the spectral response of each channel were possible, although any serious 
shifts would have been detected in the relative calibration. The sensors (all channels) were given 
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extensive absolute calibration in NASA laboratories before launch. These absolute calibrations were 
performed with the sensors and associated satellite electronics and recording systems in a simulated 
space environment. 

EVALUATION PROCEDURES 

General Procedure 

The outgoing flux densities of reflected and scattered shortwave radiation and of emitted thermal 
longwave radiation were computed as emerging from a plane and horizontal element of Earth-atmosphere 
area into the upward hemisphere. Because Nimbus 3 MRIR measured the upwelling radiances in specific 
spectral ranges only, for both types of radiation, several computational steps had to be undertaken to 
obtain from these measurements daily averages of the outgoing radiation. These steps are— 

(1 ) Computation of “total” (unfiltered) radiance in the spectral region from 0.2 to 4.0 Mm and 
4.0 to 200 Mm, respectively 

(2) Correction for the unique dependence of the measured values on the zenith and azimuth 
angle of measurements 

(3) Numerical integration over all angles to obtain the outgoing flux densities at the moment of 
measurement 

(4) Computation of the daily average of outgoing flux densities 

To perform steps (2) and (4), the geographical coordinates and the zenith and azimuth angles 
and, on the daylight side, the zenith angle of the Sun had to be computed for each single measure- 
ment. Formulas for these computations are listed in appendix A. 

To avoid errors due to shielding by heterogeneous clouds, which might not be accounted for in 
correction models used to perform step (2), all data observed with nadir angles larger than 45° were 
excluded from these evaluations. Further, all measurements of reflected solar radiation taken over 
areas where the Sun’s zenith angle f has been computed to be larger than 80° were also omitted. Under 
these conditions small errors in the satellite attitude may cause small errors in f, but the corresponding 
errors in cos are large (see eq. (6)); in addition, the reflected solar energy signal is weak in such 
situations. 

Reflected Solar Radiation 
Procedure 

Step (1), the computation of the total radiance from the. measured filtered radiance, has been 
performed simply by the assumption that the bidirectional reflectance p' computed from a measured 
radiance is representative for the entire spectral range between 0.2 and 4.0 Mm (see previous discussion): 

N'Jd', f';A, 0, 0 

p'(d\ *',r';X,0,O = — sr- 1 (6) 

S f L cos f W <P, 0 
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where 


Nj- = filtered radiance, observed over an area with geographical longitude X and latitude 0; 6' 
and 0' are zenith and azimuthal angles of observation; f ' is the Sun’s zenith angle at the 
moment of observation. 

Sy = 860.58 W m - 2 , the filtered extraterrestrial irradiance of the Sun calculated on the basis 
of Labs and Neckel’s (1968) spectral irradiance data were adjusted to a solar constant of 
S Q = 1.95 cal cm -2 min - 1 (Drummond, 1970). The spectral curve used here and those 
of other more recent investigations (Thekaekara, 1970; Drummond, 1970) coincide very 
closely. A higher accuracy in Sy can be achieved only if concurrent measurements of the 
solar constant and/or its spectral components are accomplished. 

L = d 2 /d 2 ; d and d are the mean and true Sun-to-Earth distances, respectively. 
t = the day of measurement, 
primes = values at the moment of observation. 


Steps (2) and (3), the correction for anisotropy and the integration over all angles, are merged 
into one step in the shortwave data processing by the assumption of the validity of the following 
multiplicative law, which determines the directional reflectance r (eq. (4)) at an arbitrary zenith 
angle f from the bidirectional reflectance p (eq. (6)): 


0, t) = 


Kf) 

rtf = 0 ) 


rtf = 0) 
rtf') 


rtf') 

*p(0', 0', ?') 


7rp'(0',0', r';X, 0, 0 


(7) 


The first two factors on the right-hand side of equation (7) express a normalized dependence of the 
directional reflectance of the Earth-atmosphere system on the Sun’s zenith angle. Such relations are 
shown in figure 5. The curve labeled “ocean” has been used for all observations where the reflectance 
irp' was less than 0. 10 and the equivalent blackbody temperature T b of concurrent observations of 
the radiance in the water vapor window (channel 2) was higher than 273 K. These criteria were chosen 
to select all observations over nearly or completely cloudfree ocean areas. Data from all other regions 
were processed using the curve labeled “cloud-land,” with the following exception: viz, in both hemi- 
spheres poleward of 65° latitude the curve labeled “snow” has been chosen when the reflectance up' 
was greater than 50 percent. This latter criterion has been chosen to insure that the snow relation is 
applied to bright snow or icefields only. A detailed description from published data of the derivation 
of the relations for the cloud-land and ocean curves is given in appendix B of this report. The relation 
labeled “snow” has been taken from Kondratiev (1965, table 106). In figure 5, the curve that was used 
in evaluations of Nimbus 2 data is also shown. A comparison between albedo obtained from the same 
observations but using four different methods is discussed in appendix C. 


The third factor in equation (7) performs steps (2) and (3): the correction for the anisotropic 
nature of the bidirectional reflectance of the Earth-atmosphere system and the determination of the 
directional reflectance at the time of measurement. The latter quantity, if multiplied by the incident 
solar radiation, determines the flux density of outgoing solar radiation. In applications of this third 
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0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 

cosine £ 

t= ZENITH ANGLE OF THE SUN 

figure 5.— Relations used to account for the dependence of the directional reflectance /-(f) of the Earth-atmosphere 

system on the Sun's zenith angle. 


factor, the ocean, cloud-land, and snow “surfaces” are distinct. In the former two cases a rather 
strong anisotropy of the bidirectional reflectance had to be accounted for, as it is shown in figures 6 
and 7, while diffuse reflection has been assumed over the bright polar snow and ice surfaces. The same 
angular corrections have been applied for all measurements in each of three solar zenith angle ranges; 
i.e., 0° to 35°, 35° to 60°, and 60° to 80°. 

The numerical integration over the entire daylight period of an observed area determines the flux 
density of solar radiation W' R leaving the Earth-atmosphere system over a specific area (A, 0) between 
sunrise and sunset: 


K (A, 0, t) = 




>•[?('*); A, 0, t]s' cos 


* 


( 8 ) 


where t* is an arbitrary time and t - t is the daylight period between local noon and sunset or sun- 
rise. In this integration the refraction of solar radiation has been considered at zenith angles f > 85° 
with a simple model. Without this correction, both the incoming and the reflected radiation would have 
been underdetermined by about 1 percent. S' = S Q L is the value of the solar constant S Q , modified by 
the departure from the mean of the Sun-Earth distance on day t. Daily (or 24-hr) averages of both the 
incoming and the outgoing solar radiation are 


HA X, 0, t) = 


.±r 

>2 A 


S' cos £(t*) dt* 


(9) 


1 1 




Figure 6.— Diagram of the third factor in equation (7), which is necessary to account for the anisotropic bidirectional 
reflectance in calculations of the directional reflectance from radiance measurement over land and all cloud-covered 
areas (cloud-land model). 


W' R (A, t) 

w R (k ,(p,t) = (r - t ) (10) 

12 

The ratio A = W R /H is the daily albedo (see section entitled “Definitions”), whereas the difference 
H s - W R determines the amount of solar radiation absorbed in a column of the Earth-atmosphere 
system. 
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Figure 7.— Diagram of the third factor in equation (7), which is necessary to account for the anisotropic bidirectional 
reflectance in calculations of the directional reflectance from radiance measurements over cloudfree ocean areas 
(ocean model). 


Discussion 

This data reduction procedure for albedo is based on the assumption of the validity of generalized 
reflective properties; i.e., angular dependences of the bidirectional and directional reflectances of the 
Earth-atmosphere system for all kinds of weather conditions over various surfaces except those 


13 



particular cases (polar icefields and cloudfree oceans) mentioned. These assumptions and others 
mentioned later will certainly bias the results, particularly those over various smaller scale areas. This 
is because there are no other observations yet available that describe in a complete way the reflection 
properties of the Earth-atmosphere system over various areas with different weather situations in the 
entire spectral interval between 0.2 and 4.0 pm. 

Further, the integration over the entire daylight period neglects completely any diurnal variation 
of reflection properties of the system that are due to related changes in atmosphere (e.g., cloudiness) 
and ground conditions during the course of a day. Other diurnal variations due to transient disturbances 
might be averaged out to some extent in the 1 5-day mean values used in these considerations. Over 
tropical areas, for instance, the mean lifetime of several phenomena is about 3 to 5 days (Staver et al., 
1970). Only the proper use of simultaneous observations of other satellites, such as the geostationary 
ATS spacecraft or one or two more satellites in circular orbits, would allow a rather accurate account 
of the diurnal variation; but such models have not been developed from the material available so far. 

The probable influence of the spectral sensitivity of channel 5 on the averaged reflectance has been 
discussed. Over some surfaces it might cause albedos that are slightly too low (1 to 2 percent). A com- 
parison of albedo calculated from the same data but with different models is described in appendix C. 

Outgoing Longwave Radiation 

Procedure 

Step (1), the computation of total outgoing longwave radiance, has been performed via a multiple 
least-square regression formula in which the radiances of concurrent measurements of all four infrared 
channels were related to the total outgoing longwave radiation. Such regression formulas were derived 
from radiances calculated for a set of 160 atmospheric models (10 different climatological temperature 
profiles with eight different cloud levels and having two moisture values (mostly 10 to 50 percent 
relative humidity)) for the troposphere. These models represent cloud heights and humidity con- 
ditions over almost all geographic areas. Calculations were done for zenith angles of 0°, 10°, 20°, 30°, 
and 40° with a program developed by Kunde (1967) and were checked with one by Wark et al. (1962). 

All coefficients and values of the standard error of estimate and of the multiple correlation coef- 
are listed in table 4. The general form of the regression equation for the “total” radiance N ( was as- 
sumed to be 


N t = a Q + u l N 2 + <x 2 N 2 + + « 4 A 4 + a 5 N 1 + ol 6 N 3 (11) 

where , N 2 , N 2 , and N 4 are the filtered radiances observed concurrently by channels 1 to 4. (See 
table 2.) Regression coefficients were also computed for other combinations of measured radiances 
in equation (1 1) to meet a possible failure of one of the four channels. Most emphasis in equation (11) 
was placed on measurements of channel 2, because its measurements were expected to be more highly 
correlated to variations of the outgoing total radiance than those of other channels. 

The contributions of the channel 2 and 4 terms in equation (11) determine nearly 96 percent of 
the variable part of the total radiance, while the rest is contributed by measurements of the other two 
infrared channels. This treatment of the measurements of all four infrared channels as four independent 
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pieces of information is not strictly correct because under many meteorological conditions data from 
a certain channel may contain much information from overlapping areas of the atmosphere covered 
by another channel. However, together they give the best possible parametrization of infrared radia- 
tion to space. 

Steps (2) and (3), the correction for angular dependence and the computation of the outgoing 
flux density at the moment of measurement, were performed in one step following the procedure de- 
veloped by Wark et al. (1962). When the emitted longwave radiance is assumed to be independent of 
the azimuth angle \jj, equation (2) for the outgoing longwave radiation W L assumes the form 

rn/2 

W L (X, (j), t) = 2n I N { (6, X, 0, t) cos 0 sin 0 d9 (12) 

Jo 


where N is the total radiance obtained by use of equation (11) from “filtered” radiance measurements 
over an area having the geographical coordinates X and 0. Assuming a generalized “limb-darkening 
function” /(0), which describes the change of radiance with the zenith angle of measurement, equation 
(12) can be simplified to 


where 


and 


W L (X,<l>,t) = N t (d= 0;X, 0, t)Y 


Y-2tt 



f(6) cos 0 sin 0 dd 


N t W) N s ^(d) 

f(6 )= * 

N t {8 = 0) )V 5 (0 = 0) 


(13) 

(14) 


(15) 


Here 0 is the zenith angle of measurement. Limb-darkening functions /(0) were derived from radiance 
measurements of Nimbus 2 (Nimbus Project Staff, 1966) taken in the spectral region between 5 and 
30 jum over the entire globe in 12 days during the period from May 16 to July 29, 1966. To observe 
different climatological conditions, data from five different areas were studied; their coordinates and 
some climatological characteristics of their Earth-atmosphere system are listed in table 5. Data of 
areas 1 to 4 were grouped into two samples of six successive days to test for a seasonal trend of the 
limb darkening. In addition, only measurements made at nadir angles of less than 54° were used to 
avoid “space contamination.” These data were then fitted by a least-square method to the following 
equation: 

A 5 - 3O (0) = /V 5 . 30 (0 = OK 1 +j3,0+j3 2 0 2 +0 3 0 3 ) (16) 


The Nimbus 2 average values of N 5 3Q (0 = 0), the Y integrals, the standard errors of estimates, and 
standard deviations of N 5 3O (0) are summarized in table 6. 

The plots in figure 8 of the limb-darkening function /(0) at zenith angles of less than 45° (the 
limit set for radiation budget use of the measurements) fall into two classes. Thus in the final evalua- 
tion procedure only two different relationships for /(0) (eq. (15)) have been used to discriminate 
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Table 5.— Geographical Areas of Nimbus 2 Data Samples Used To Derive Limb-Darkening Functions f{Q) 


Number 

Area 

Location 

Notes 

1 

Arctic 

70° N < 0 < 90° N, 
0° < X < 360° 

Low surface temperature (—270 K), 
dry, relatively high stratosphere 
temperature (~235 K) 

2 

Desert 

10° N < 0 < 35° N, 
310° W< A <360° 

High surface temperature (~3 1 0 K), 
dry, low stratosphere temperature 
(-205 K) 

3 

Tropical and 

30° S < 0 < 35° N, 

High surface temperature (—295 K), 


subtropical 

ocean 

110° W < X < 200° W 

dry and moist, high and cold 
clouds 

4 

Southern 

60° S < 0 < 30° S, 

Low surface temperature (—280 to 


midlatitudes 

0° < X < 360° 

285 K), moist, very cloudy, cold 
stratosphere (—222 K) 

5 

Antarctic 

90° S <0 < 70° S, 
0° < X < 360° 

Very low surface temperature 
(—210 to 230 K), dry, cold 
stratosphere (—215 K) 


Table 6. -Infrared Limb Darkening of the Earth- 
Atmosphere System 


Area 

N S-3q(P = 0), 
W m~ 2 sr 1 

Y, sr 

— 
Standard error 
of estimate, 
Wm' 2 sr -1 

1 

48.9 

3.05 

3.7 

la 

46.4 

3.05 

3.8 

2 

69.7 

2.96 

10.4 

2a 

68.0 

2.80 

11.2 

3 

61.3 

2.93 

8.3 

3a 

61.2 

2.89 

8.0 

4 

48.7 

2.93 

8.0 

4a 

49.3 

2.92 

8.0 

5 

29.3 

2.95 

6.4 


between polar (poleward of 70° latitude) and remaining areas. The two short-dashed curves for the two 
desert samples (2, 2a) in figure 8 deviate from each other considerably. The only explanation that can 
be offered is that the relatively small area encompassed by the Sahara and Arabian deserts severely re- 
stricted the sample population compared to the other regions. In the analysis the mean tropical and 
midlatitude limb-darkening curve (from 3, 3a, 4, and 4a) was used for low-latitude and midlatitude 
deserts. 

Step (4), the calculation of long-term daily averages, was performed finally with 1 5-day averages 
of the outgoing longwave radiation obtained in each grid field from day and from night measurements, 
where each value has been weighted according to the length of the daylight and nighttime periods, re- 
spectively. If a grid field lacked one value, either day or night, no daily average has been obtained, 
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Figure 8.— Limb darkening of the Earth-atmosphere system in the spectral region between 
5.0 and 30.0 pm as obtained from Nimbus 2 measurements over various areas listed in 
table 5. 


thus, no value of the net radiation Q could be determined for that area. In the final map analysis, such 
gaps were given special notation. 

Discussion 

In the approach used here to determine the limb-darkening functions, it is assumed that the Earth- 
atmosphere system is uniform in its horizontal temperature distribution and is a plane over each area 
of consideration, because this procedure contains no automatic method for distinguishing different 
cloud levels or cloud-covered from cloudfree areas. The standard error of estimate in table 6 is very 
high because of the scattering of available data. Thus, these limb-darkening functions can be con- 
sidered to represent only average conditions valuable for calculations of the outgoing longwave 
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radiation over longer periods. A more specific treatment requires more sophisticated observations such 
as might be provided in a few years from the Earth radiation budget experiment to be flown on Nimbus 
F. 


This entire evaluation procedure encounters several independent error sources of which step (1), 
the calculation of the total radiance N from four independent measurements, might be the most im- 
portant one. Two error sources are possible in this case: random noise in each channel and systematic 
trends in one or more channels. 

The noise in all four infrared channels originates primarily by a digital number error of 1.0 in the 
MRIR-A/D converter on board Nimbus 3. It causes errors in the equivalent blackbody temperatures 
T b of more than 1° only if T b < 230 K (except channel 2: T b <210 K). Its values are estimated to be 
0.007, 0.07, 0.05, and 0.04 W m~ 2 sr -1 for channels 1 through 4, respectively (McCulloch 1 ). These 
noise values and multiples of them were introduced into equation (1 1) as a random process. They in- 
fluenced the standard error of estimate of the regression procedure by less than 1 percent. Indeed, 
when the noise level has been increased over actual estimates by a factor of 40, the standard error of 
estimate increased only by a factor of 10. 

Errors may also arise if one or more of the four channels are degrading in their total and spectral 
sensitivity. Degradations in the total sensitivity were not observed and could be checked easily with 
the onboard calibration system. However, there was no means to check the spectral sensitivity. 

Another error source, which has not been considered yet, may be caused by systematic errors in 
the model calculations. A complete modeling of infrared radiation under natural conditions is very 
difficult. Thus, results of the present study potentially contain more serious errors than the evaluation 
of Nimbus 2 data (Raschke and Bandeen, 1970), because here four different and narrow spectral in- 
tervals had to be considered. For Nimbus 2 only measurements of a very broad channel (5 to 30 pm) 
had to be converted into total radiances. If independent and concurrent measurements of the upwell- 
ing total radiance from the same field of view were available for correlation with the other four 
measurements, then one should obtain a somewhat poorer correlation than obtained by calculation, 
but a more accurate parametrization of natural conditions. Finally, the results are also not satisfactorily 
corrected for diurnal variations of temperature and cloudiness, although over each area, measure- 
ments during daylight and nighttime conditions were available (primarily local noon and midnight). 
Figure 9 demonstrates the resulting error in radiation balance caused by errors in the albedo and the 
outgoing longwave radiation. 

RESULTS 

Global and Hemispherical Radiation Budget 

All global and hemispherical averages of the radiation balance Q and its three components as ob- 
tained from all Nimbus 3 MRIR measurements within semimonthly periods are listed in table 7. Annual 
averages of the planetary albedo and the outgoing longwave radiation for the globe and hemispheres 
were obtained from a graphical presentation of their respective semimonthly averages; the annual 


1 A. W. McCulloch, private communication, 1969. 
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dH s = °-° A = ALBEDO 

H s =0.2 cal cm' 2 min -1 

H s = 0 .5 cal cm' 2 min-' 


Figure 9.— Error diagram for estimates of the absolute error in the radiation balance Q for two 
known irradiances of incoming solar radiation H . Error caused by errors in the albedo A and in 
the radiant exitance of outgoing longwave radiation W L . 


averages of the incoming solar radiation were computed directly for the solar constant. All numbers 
in table 7 are written with an accuracy of better than 1 to 2 percent, although it is believed in this 
study that the third decimal should be considered to be not very accurate. 

These results show that during a year the incoming solar radiation of 0.488 cal cm - 2 min - 1 is 
almost completely balanced by a global albedo of 28.4 percent and a value of outgoing longwave 
radiation of 0.345 cal cm - 2 min - 1 . This latter value corresponds to a mean blackbody temperature of 
- 18° C or 255 K. These global annual averages confirm in their magnitude older results found by Vonder 
Haar (1968) and Vonder Haar and Suomi (1971) from measurements of several other satellites that 
carried hemispherical and flatplate radiometers. 

All previous investigations of the radiation balance with climatological data, however, resulted in 
a much lower global emission, corresponding to a mean blackbody temperature of ~250 K and a 
higher global albedo of more than 33 percent. A comparison of respective zonal averages will reveal 
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Table 7.— Global and Hemispherical Radiation Budget of the Earth- Atmosphere System 


Date 

Solar radiation, 
cal cm - 2 min - 1 

Albedo 

Outgoing 
longwave 
radiation, 
cal cm -2 min -1 

Radiation 
balance, 
cal cm - 2 min - 1 

Incoming 

Absorbed 

Apr. 16 to 30, 1969: 






N 

0.586 

0.410 

0.300 

0.349 

0.061 

S 

.379 

.277 

.270 

.342 

-.065 

G 

.483 

.344 

.288 

.346 

-.002 

May 1 to 15, 1969: 






N 

.618 

.431 

.302 

.351 

.080 

S 

.334 

.244 

.269 

.343 

-.099 

G 

.476 

.337 

.291 

.347 

-.010 

May 16 to 31, 1969: 






N 

.641 

.440 

.314 

.355 

.085 

S 

.310 

.228 

.264 

.346 

-.118 

G 

.476 

.334 

.298 

.350 

-.016 

June 1 to 15, 1969: 






N 

.656 

. .454 

.308 

.359 

.095 

S 

.290 

.213 

.267 

.343 

-.130 

G 

.473 

.333 

.296 

.351 

-.018 

June 16 to 30, 1969: 






N 

.659 

.462 

.299 

. .362 

.100 

S 

.283 

.206 

.271 

.342 

-.136 

G 

.471 

.334 

.291 

.352 

-.018 

July 1 to 15, 1969: 






N 

.652 

.459 

.296 

.365 

.094 

S 

.291 

.213 

.267 

.344 

-.131 

G 

.471 

.336 

.286 

.355 

-.019 

July 16 to 31, 1969: 






N 

.633 

.449 

.290 

.363 

.086 

S 

.311 

.230 

.261 

.345 

-.115 

G 

.472 

.339 

.281 

.354 

-.015 

Aug. 1 to 15, 1969: 






N 

.609 

.434 

.287 

.365 

.069 

S 

.340 

.248 

.271 

.341 

-.093 

G 

.474 

.341 

.280 

.353 

-.012 

Oct. 3 to 17, 1969: 






N 

.443 

.323 

.270 

.351 

-.028 

S 

.541 

.384 

.291 

.344 

.040 

G 

.492 

.353 

.282 

.347 

.006 

Jan. 21 to Feb. 3, 1970: 






N 

.343 

.249 

.273 

.332 

-.083 

S 

.663 

.473 

.287 

.342 

.131 

G 

.501 

.361 

.283 

.337 

.024 

Annual: 






N 

.483 

.344 

.287 

.346 

-.002 

S 

.492 

.354 

.280 

.344 

.010 

G 

.488 

.349 

.284 

.345 

.004 


N = northern hemisphere, S = southern hemisphere, G = globe. 
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Figure 10.— Annual zonal averages of albedo (percent) and outgoing longwave radia- 
tion obtained from Nimbus 3 (dashed line) and earlier satellite (solid line: Vonder 
Haar and Suomi, 1971) and from calculations with climatological data (open circles: 
London, 1957). 


in greater detail these differences and provide some basis for discussion of possible error sources. (See 
fig. 10). 

The northern hemisphere receives during a year per area and time unit about 2 percent less solar 
radiation than the southern hemisphere, because of the characteristics of the Earth’s orbit around the 
Sun. Its albedo, probably due to the brighter land surfaces and larger snow-covered areas, is somewhat 
higher than that of the southern hemisphere over which, however, the cloudiness might be somewhat 
greater (Van Loon, 1970); but there the mean surface albedo is very low due to a preponderance of 
ocean areas. Thus, the Earth-atmosphere system in the northern hemisphere absorbs somewhat ( — 2.5 
percent) less radiative energy than that over the southern hemisphere. It was found from these data 
that its longwave emission is slightly higher than that of the southern hemisphere. 

Thus, provided these results are correct in their relative magnitudes during this particular period, 
some horizontal circulation processes must have transported energy from the southern to the northern 
hemisphere, which was determined from these measurements to be a deficit area. The absolute magni- 
tude of the radiative balance over both hemispheres should be almost equal to obtain global balance 
Q = 0. The small value of Q = 0.004 obtained in these investigations fulfills almost completely this re- 
quirement. However, in an objective discussion of these results in table 7, their limitations should be 
kept always in mind: 

(1 ) The results were obtained from a still very incomplete data set. 

(2) The evaluation methods might bias the results in a yet indeterminable direction. 

(3) Data were collected during only one single annual period. 

(4) An assumed value of the solar constant was used. 
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Small seasonal variations of the order of 2 to 4 percent of the annual averages can be observed in 
the results for the global albedo and the outgoing longwave radiation. The albedo is slightly higher 
during the period from April to June than in the remaining part of the year, possibly due to the 
brighter land surfaces and their snow-covered parts of the northern hemisphere which are missing over 
the southern hemisphere. The outgoing longwave radiation over the northern hemisphere shows a 
small seasonal variation, probably related to the heating and cooling of its continents, while over the 
southern hemisphere it is almost constant throughout the entire year because of the preponderance of 
ocean surfaces. The trends of global averages of the absorbed and outgoing radiation during the period 
May to July can also be observed in results from Nimbus 2 measurements (Raschke and Bandeen, 1970). 
However, the values in this report are 1 to 2 percent higher and lower, respectively, than those earlier 
results, primarily because of the different empirical models used in the computations. The global radia- 
tion balance Q shows a slight deficit during the northern hemispheric spring and summer and a gain of 
radiative energy in both other seasons, due to the previously mentioned seasonal changes and also the 
change in the distance between the Sun and Earth. Variations of Q of a similar phase, but higher mag- 
nitude, were estimated by Simpson (1929) from observations available at that time. The entire globe 
is close to radiative equilibrium during all months within the assumed accuracy of the measurements, 
although the incident solar radiation changes its value by 7 percent throughout the year due to the 
varying distance between Sun and Earth. 

Zonal Radiation Budget 

In contrast to global hemispherical values, the zonal averages of the radiation balance and its 
components, shown in time-latitude diagrams in figures 11 to 14, have a pronounced seasonal variation 
in almost all latitudes. Many details can be observed in the results from all measurements during the 
period from April 16 to August 15, 1969, while over the remaining and larger part of the year some 
could be obtained by various and not necessarily valid interpolations. 

In tropical and subtropical latitudes of both hemispheres the variations of the albedo and out- 
going longwave radiation follow closely the mean global movement of the intertropical convergence 
zone (ITCZ) and associated changes of the cloudiness in the subtropics. The mean global ITCZ can be 
identified with a mean albedo between about 25 and 30 percent and an emission of less than 0.35 
cal cm" 2 min - 1 . In the northern midlatitudes the albedo decreases rapidly from winter to summer, 
because of snow melt; also the continents heat considerably, having the highest emission in July and 
August. However, only small corresponding changes occur over the southern hemisphere at the same 
latitudes. 

Over the Arctic regions, poleward of 70° N, the albedo is on the average not higher than 70 per- 
cent, while from Nimbus 2 measurements even higher values were found (Raschke and Bandeen, 

1970) because another relation for integration over the entire daylight period was used. (See subsection 
entitled “Reflected Solar Radiation” and fig. 5.) In winter (January) the emission of 0.23 cal cm -2 
min -1 corresponds to temperatures as low as 230 K (-43° C), while in summer (July) these equivalent 
blackbody temperatures are 20° to 30° higher. During the polar night over the Antarctic Continent, 
the emission temperatures fall to values as low as 218 K (-55° C) (in July), while in summer these tem- 
peratures are almost 20° higher. 
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Figure 13.— Time-latitude diagram of the 
measurements in 1969 and 1970. Dashe 
wave radiation. 
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Figure 14.— Time-latitude diagram of the radiation balance (calories per square centimeter per minute) determined from Nimbus 3 measure- 
ments in 1969 and 1970. Dashed portions of isolines were determined by interpolation. Shaded portions are regions of negative net 
radiation. 




For the reader’s convenience, in appendix E a tabulation of the Stefan-Boltzmann law is presented, 
providing a simple method for estimation of an effective emission temperature from values of the out- 
going longwave radiation. 

The albedo of the Antarctic ice, as found from only two semimonthly periods of measurements, 
is not higher than 60 to 75 percent. In the central Antarctic it decreases from 70 to 60 percent between 
October and January, probably due more to an insufficient correction for the change of the mean insola- 
tion angle than to surface changes, such as snow and ice cover and their reflection properties. All snow 
and ice albedos appear to be too low, but they are values considered to represent the entire spectral 
interval between 0.2 and 4.0 /xm. According to the earlier discussion, they might have been under- 
estimated by a factor of about 1.03. 

The seasonal change of the radiation balance follows closely, but with some phase delay, the 
Sun’s declination, which is shown in figure 14. Areas of major surplus of radiative energy are in each 
corresponding season the subtropics, while the major deficit is found poleward of about 60° latitude. 
Over the southern hemisphere a reversal of the gradient of the radiation balance is found at all seasons 
because of the very low temperatures above the high plateaus of the Antarctic Continent. There, even 
in summer the balance shows a deficit of more than 0.04 cal cm - 2 min - 1 , while in June and partly in 
July 1969 the Arctic ice shield obtains slightly more radiation than it loses to space. This was also 
found from earlier satellite measurements (Vonder Haar, 1968). 

Geographical Distributions of the Radiation Balance During the "High" Seasons 

Maps of the geographical distribution of the radiation balance, the albedo, the outgoing longwave 
radiation, and the absorbed solar radiation were produced from measurements of each semimonthly 
period and analyzed on Mercator (between 40° N and 40° S) and polar stereographic (poleward of 40° 
latitude) projections. The spatial resolution in these ranged from about 500 by 500 km 2 (between 20° 

N and 20° S) to better than 250 by 250 km 2 (poleward of about 40°) at each gridpoint. 

Some of these maps were transferred manually and by a computer program into Aitoff projec- 
tions, which allow a convenient but less resolved view on the patterns of each quantity over the entire 
Earth. Maps for the “high” seasons (July 16 to 31, 1969, and January 21 to February 3, 1970) are dis- 
cussed in this section. Those of two other periods (May 1 to 15 and October 3 to 17, 1969) that may 
represent the spring and fall seasons are shown in appendix D. 

These geographical distributions of the albedo (figs. 1 5 and 1 6), the absorbed solar radiation (figs. 
17 and 18), the outgoing longwave radiation (figs. 19 and 20), and the radiation balance (figs. 21 and 
22) reveal clearly the major areas of gain and deficit of radiative energy. The overwhelming latitudinal 
gradient in poleward direction is accompanied by longitudinal gradients caused primarily by the land- 
sea distribution. These weaker gradients may also, as does the poleward gradient, cause mechanisms of 
energy exchange. 

One outstanding feature in this pattern is the deficit found in all seasons over the bright (albedo 
> 35 percent) and hot W L > 0.39 cal cm - 2 min - 1 ) desert regions of North Africa and the Arabian 
Peninsula. The ocean areas at the same latitude are the major areas of heat storage in the spring, sum- 
mer, and fall. 
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Figure 15.— Albedo (percent) of the Earth-atmosphere system during the period July 16 to 31, 1969. Values are daily averages. 




NO ILLUMINATION 



Figure 16.-Albedo (percent) of the Earth-atmosphere system during the period January 21 to February 3, 1970. Values are daily averages. 
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Figure 17.— Solar radiation (calories per square centimeter per minute) absorbed in the Earth-atmosphere system during the period July 16 to 31, 1969. 

Values are daily averages. 
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Figure 18. -Solar radiation (calories per square centimeter per minute) absorbed in the Earth-atmosphere system during the period January 21 to 

February 3, 1970. Values are daily averages. 
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Figure 19.— Outgoing longwave radiation (calories per square centimeter per minute) emitted from the Earth-atmosphere system to space during the 

period July 16 to 21, 1969. Values are daily averages. 


NO DATA 
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Figure 20.— Outgoing longwave radiation (calories per square centimeter per minute) emitted from the Earth-atmosphere system to space during the 

period January 21 to February 3, 1970. Values are daily averages. 
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Values are daily averages. 
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Figure 22-Radiation balance (calories per square centimeter per minute) of the Earth-atmosphere system during the period January 21 to February 3, 

1970. Values are daily averages. 




Elsewhere in midlatitudes the patterns in the albedo and longwave emission maps are related 
mostly to dynamical processes (bright and mostly cool clouds; or dark, cloudless, and warm ocean 
areas) and the seasonal temperature patterns. In January the albedo of the Antarctic icefields, as in 
July over Greenland, does not exceed values of 65 percent, which is considerably lower than it has 
been reported in the literature from ground observations (80 percent and more; Hoinkes, 1968) and 
from Nimbus 2 observation (70 to 80 percent; Raschke, 1968). 

These differences between albedo values at the ground and at the top of the atmosphere are 
mainly due to the fact that the atmosphere absorbs completely the incident near-infrared radiation 
contained in the Sun’s extraterrestrial spectrum. For ground-based measurements of both the reflected 
and incoming solar radiation, this component is very small. Thus, these measurements result in albedo 
values that apply primarily to the visible and nearest infrared, while the albedo at the top of the atmos- 
phere is an average over a much wider spectral range. 

CONCLUSIONS 

The principal task of this report has been to discuss the evaluation method that was used to com- 
pute, from Nimbus 3 measurements of reflected shortwave (0.2 to 4.8 pm) and emitted longwave (four 
narrowband channels between 6 and 24 pm) radiation, the radiation balance of the Earth-atmosphere 
system, and to discuss its error sources. As examples of results, maps are included presenting the geo- 
graphical distributions of the albedo, the absorbed solar radiation, the outgoing longwave radiation, 
and the radiation balance results that were obtained from 1 5-day periods of measurements during the 
four seasonal periods (May 1 to 15, July 16 to 31, and October 3 to 17, 1969, and January 21 to 
February 3, 1970). In addition, global, hemispherical, and zonal averages of the radiation budget param- 
eters are presented for all Nimbus 3 1 5-day measurement periods. 

In summary, annual global averages of the albedo of 28.4 percent and of the corresponding out- 
going longwave radiation of 0.345 cal cm - 2 min - 1 were obtained. The annual incoming solar radiation, 
together with these values, balances to within less than 1 percent (+0.004 cal cm -2 min - 1 ) of the 
annual global solar irradiance (which has been computed for a solar constant of S Q = 1 .95 cal cm - 2 
min - 1 ). These results confirm those of earlier investigations of the Earth’s radiation budget with satel- 
lite measurements, which indicated a darker and warmer planet Earth than previously found from 
numerical studies with climatological data. 

Principal error sources in the evaluation technique occur in the various steps explained in the sec- 
tion of this report entitled “Available Data” and in appendix B. Primarily the calculation of the total 
(4.0 to 200.0 pm) radiance from measured filtered radiances of outgoing longwave radiation and the 
use of gross-empirical models for the calculation of daily averages of outgoing radiant flux densities of 
solar radiation from radiance measurements can bias the results. Unfortunately, no conclusive com- 
parisons of simultaneous observations with radiometer and albedometer sondes were available to check 
these results. However, the agreement with results from satellite observations during previous years is 
good. It may be concluded from this study that the absolute accuracy of these results may be only 
within 5 percent of the albedo and the outgoing longwave radiation. 

Further investigations should include simultaneous observations of the flux density of incident 
solar radiation, because its value has been assumed constant. Much more data are needed to establish 
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better and more representative empirical models on the reflection properties of the Earth-atmosphere 
system. Consequently, extensive checks of results are needed with in-situ observations from either 
balloon or airplane experiments. 

These results, like those of earlier satellite measurements, cannot be used for studies of long-term 
climatic changes of natural or manmade origin because their evaluation was based on too many assump- 
tions. Once better models on the reflection properties of the Earth-atmosphere system become avail- 
able, a reevaluation of these data would appear to be advisable. More representative results on the time 
and space spectrum of the Earth’s radiation budget can be obtained from simultaneous measurements 
with scanning and integrating (flux) radiometers on board several satellites whose orbital characteristics 
meet the requirements of accurate data sampling. 
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Appendix A 

COMPUTATION OF ANGLES ?, Q, AND 0 


SOLAR ZENITH ANGLE f 

f = cos -1 [cos 4> cos 8 cos (X G - X) + sin 0 sin 5] (A-l) 

where X and 0 are the geographic longitude and latitude, respectively, of an observed surface element 
(X is counted westward; 0 is positive north and negative south of the equator), X G is the Greenwich 
hour angle, and 5 is the declination of the Sun. 


ZENITH ANGLE OF MEASUREMENT 9 

6 = sin" 1 1 — — — sin 

= sin" 1 ( K sin a) (A-2) 

where a is the nadir angle of measurement, which can be obtained from the data tape. The Earth’s 
mean radius/? is assumed to be 6371 km and the satellite height //to be 1120 km; therefore, 


K = 


R + H 
R 


1.1758 


AZIMUTH ANGLE OF MEASUREMENT 0 

The azimuth angle of measurement 0 is determined with respect to the ray of incident solar 
radiation (fig. A-l) by 


where 


0 = cos 1 


cos f cos 6 - cos T 
sin f sin 6 


(A-3) 


[ K cos 0 y cos (X s - \ G ) - cos 0 cos (X - X c )] cos 5 + (K sin 0 s - sin 0) sin 6 

cos p = ■ — — - ' (A-4) 

\K 2 + 1-2 K [cos <t> s cos 0 cos (X - Xp + sin 0 f sin 0] 

where X s and 0 S are the geographic longitude and latitude, respectively, of the subsatellite point. 
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ZENITH 



Figure A-1.— Geometrical configuration of Earth and Sun. 


LENGTH OF A DAY AX 


Sunrise and sunset are defined as the instants when the upper edge of the disk of the Sun is on 
the horizon at normal refraction (table IX in Baur, 1953). Assuming that a semidiameter of the Sun 
subtends a 16-min angle with a constant 34-min refraction (List, 1963), then the length of a day AX* 
(in degrees of arc on the celestial sphere) is computed as follows: 


AX* = 2AX = 2 cos" 1 


0.0143 - sin <p sin 5 
cos (j> cos 5 


= 2 cos 1 (arc) 


(A-5) 


where AX is 1 80° ( 1 2 hr) if (arc) < - 1 . 
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Appendix B 


DERIVATION OF INTEGRATION MODELS FOR THE CALCULATION OF 
THE OUTGOING SHORTWAVE RADIATION 


DATA AND GENERAL OUTLINE 

After the evaluation of Nimbus 2 radiation measurements (Raschke, 1968) was made, more 
data (both measurements and calculations) were published and thus became available for use in a 
systematic consideration of the anisotropic reflection properties of the Earth-atmosphere system. 

These data allowed a distinction to be made between the reflection properties of cloudfree ocean 
areas and of cloud-covered and land areas. They still cannot be considered to satisfy completely the 
requirements necessary for an accurate calculation of the flux density of outgoing shortwave radiation. 
The procedure follows strictly that previously used for evaluations of Nimbus 2 data. 

The sources of data used in these derivations of numerical integration models are summarized 
in table B-l. 

This material was very heterogeneous with respect to its origin, spectral and angular range of 
observation, and area of measurements. Thus, another multiplicative law (eq. (B-l)) had to be 
assumed to obtain a mean curve for the dependence of the directional reflectance r(f) on the Sun’s 
zenith angle: 

/-(?) _ 7rp( 0, 0, n /-a) 7TP(0, 0, 0) 

z(? = 0) ttp( 0, 0, 0) 7rp(0, 0, ?) rtf = 0) ( ' 

The three arguments of p are 6, <//, and f. In the first factor of equation (B-l) the change of the 
bidirectional reflectance toward the zenith (or observed at the nadir) with the zenith angle of the 
Sun f is considered. The second factor and the third factor relate the directional reflectance r(f) to 
the bidirectional reflectance toward the zenith. The use of equation (B-l) implies the physical 
assumption of a unique relation as expressed by these three factors for all natural conditions 
(different surfaces and/or various clouds at various altitudes) in the Earth-atmosphere system. 

MODELS FOR CLOUDY CONDITIONS AND FOR LAND SURFACES (CLOUD-LAND) 

This specific model has been applied to all measurements except those taken over cloudfree 
oceans (channel 2: T b > 273 K and 7 rp' < 10 percent) and solid snow or ice masses <p > 65° latitude 
and 7rp' > 50 percent). The relation 




p(Q, o, D 
p( o, o, 0) 


(B-2) 
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Table B-l.— Data for Derivation of Empirical Reflection Models 


Reference 

Type of measurement or analysis 

Surface 

Bartman 

MRIR on balloon (0.2 to 4.0 pm) 

Snow 

(1967) 



Salonionson 

MRIR in small airplane (0.2 to 4.0 

Snow stratus 

(1968) 

pm) 


Cherrix and Sparkman 

MRIR in Convair 990 (0.2 to 4.0 pm) 

Ocean surface, stratus, and 

(1967) 


broken clouds 

Griggs and Marggraf 

Airplane with albedometer 

Clouds and ocean surface 

(1968) 



Brennan (1969) 

MRIR in Convair 990 (0.2 to 4.0 pm); 
some full rosettes, but mostly 
measurements in the principal 
plane and at \p = 90° 

Ocean 

Brennan and Bandeen 

MRIR in Convair 990 (0.2 to 4.0 pm); 

Surface 

(1970) 

some full rosettes, but mostly 
measurements in the principal 
plane and at i p = 90° 


Ruffetal. (1968) 

TIROS 4 MRIR (0.55 to 0.75 pm) 
statistical analysis 

Clouds 

Arking 3 

TIROS 4 MRIR (0.55 to 0.75 pm) 
statistical analysis 

All types of surfaces 

Raschke and Bandeen 

Photomultiplier in ATS 1 

Cloudfree oceans 

(1969) 

(0.45 to 0.65 pm) 


Korb and Moller 

Calculations for plane clouds 

Clouds 

(1962) 



Plass and Kattawar 

Monte Carlo calculations at various 

Smooth ocean 

(1968) 

wavelengths 



a A. Aiking, GSFC, private communication, 1967. 


has been derived from the statistical analyses of various authors, as shown in figure B-l. These results 
show considerable disagreement at low Sun (J > 70°) and high Sun (f < 40°). For the derivation of a 
mean curve, the values of all three investigations were averaged for f < 70°. At a very low Sun two of 
the investigations shown in figure B-l tend to show values of F x <1, thus it has been assumed to 
be 0.9. 

The second and third factors of equation (B-l) were obtained from all data summarized in 
figure B-2, whose ordinate is 

F m = '(» n P(®> 0, 0) 

2 5) = 0) 77 P (0, 0, f ) 

= r«)/7rp(0, 0, g) 3 , 

■ K? = 0)/ttp(0, 0, 0) 
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P (0, 0,0 



ZENITH ANGLE OF THE SUN 


Figure B-1 .-The function F 1 (f) = p( 0, 0, f)/p( 0, 0, 0) versus the Sun's zenith angle. 



10° 20° 30° 40° 50° 60° 70° 80° 90° 

ZENITH ANGLE OF THE SUN 


Figure B-2.-The function F 2 (f) = [r{f)//-(f=0)l [7rp(0, 0, 0)/vp(0, 0, versus the Sun's zenith angle for 

cloud- and snow-covered areas. 
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The mean curve in this figure follows almost completely the results by Ruff et al. (1968). Values of 
F 2 as shown in figure B-2 are required to construct mean diagrams r($)/irp(d, \p, f) for the calculation 
of the directional reflectance from an “observed” bidirectional reflectance as described in the section 
entitled “Reflected Solar Radiation.” Curves obtained from observations and used to construct these 
diagrams (fig. 6) are shown in figures B-3, B-4, and B-5. The curves in figure B-3 show that all 
measurements taken over areas at high Sun (0 < f < 35°) need almost no correction for the anisotropic 
angular reflection. Finally the curve r(f)/r(f = 0) that describes the change of the directional reflectance 
with the Sun’s zenith angle f is obtained by multiplication of the two curves derived in figures B-l and 
B-2. This curve is shown in figure B-6 (see also fig. 5) and compared with some results that Korb and 
Moller (1962) obtained by model calculations for thick cumulus and stratus layers. These calculated 
results show a less steep slope, while others by Plass and Kattawar (1968) follow this curve very 
closely. They, however, were determined for one specific wavelength only, which is described here 
by an optical thickness r = 10, and for a very-low-reflecting surface ( A = 20 percent). 

MODELS FOR CLOUDFREE OCEAN AREAS (OCEAN) 

Assuming that the optical properties of the ocean-atmosphere system did not deviate from each 
other at each measurement, a mean curve of r(J) could be obtained by interpolation and hand 
smoothing through a collection of data by Griggs and Marggraf (1968) and by Brennan (1969). This 
curve is shown in figure B-7. Brennan also published data of two full rosettes taken at two different 
angles of illumination by the Sun (29.5° < f <35.5°; and 44.6° < f < 50.5°). They were used to 
construct the diagram of r(£)/np(d, \p, f) shown in figure 7. 

The specular reflection of direct solar radiation causes most of the anisotropy at 6 < 50° when 
the Sun is higher than f = 50°. (See also Raschke, 1971 .) Thus, there is almost no anisotropy 
correction needed in this range of 6, if f < 60°. Most of Brennan’s measurements, unfortunately, 
were obtained in the principal plane only. The directional reflectance r(f) has been calculated from 
them assuming that the bidirectional reflectances observed in the forward direction (\jj = 0°) are 
representative for an angular range of A\jj = 30° centered at \p = 0°, because the Sun’s glint on the 
water surface covers only a small angular range. The backward measurements (i p - 180°) are 
representative for all other angles of \p. This method is very debatable, but it seemed to be the only 
way to estimate the directional reflectance of the oceans from these measurements. The curve is 
shown in figure B-8. 

DISCUSSION 

As stated earlier, the available set of data is very heterogeneous and incomplete. Under these 
conditions, one must expect systematic errors whose magnitude and direction (too low or too high 
albedo) cannot be estimated unless there are proper measurements available. Such measurements, 
however, are not available, and thus a complete error analysis is not possible. In particular, data on 
bidirectional reflectance characteristics of snow and ice were so very sparse that a nearly diffuse 
pattern was used in the data reduction (i.e., r(^)/7rp(0, i (/, » 1.0). 
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Figure B-3.— The ratio r(f)/7rp(0, f) obtained from various observations at 0° < f < 35' 
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Figure B-4.-The ratio r(i;)/-np(6, i//, f) obtained from various observations at 35° < f < 60°. 
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Figure 8-5. -The ratio r{$)hip{6, \jj, f) obtained from various observations at 60° < f < 80°. 



2.0 


0 oKORB AND MOLLER (1962), cu (1000 m) 

» ® KORB AND MOLLER (1962), sf (500 m) 



Figure B-6.— The ratio r(f)/r(f = 0); cloud-land model. 



Figure B-7.— Directional reflectances /■($■) as measured over the cloudfree ocean. 
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Appendix C 


COMPARISON OF THE DIFFERENT MODELS USED TO ACCOUNT 
FOR THE INCREASING DIRECTIONAL REFLECTANCE r(f) 
WITH INCREASING ZENITH ANGLE OF THE SUN 


The snow, cloud-land, and ocean models show (fig. 5) an increasingly steep slope. A qualitative 
explanation of this behavior is that the models (except for the Nimbus 2 model) originate from 
measurements over surfaces with very high, mean, and very low albedo. Qualitatively equal results 
can be found from multiple scattering calculations of the field of shortwave radiation in various 
atmospheric models with ground surfaces of different albedos because the flux density of shortwave 
radiation at the top of the atmospheric model is increasingly affected by multiple scattering in the 
atmosphere with decreasing reflectance of the surface. 

Each of the models in figure 5 might result in a calculated value for albedo that is too low or too 
high if used with the wrong data; therefore, in this appendix a comparison of albedo values is discussed. 
They were obtained from a data sample (orbit 861 primarily) applied to different models. The 
Nimbus 2 model (fig. 5) was originally derived from various information available in 1967 when the 
Nimbus 2 data evaluations were being made. Thus, it represents a weighted mean of all results 
available at that time. 

Figure C-l shows that the albedo over the Greenland icepack and over the Antarctic would be 
overestimated if the cloud-land instead of the snow model were used, especially at the highest magni- 
tudes. Over the brightest areas, having an albedo of about 70 percent according to the snow model, 
one would have obtained albedos of more than 80 percent if the cloud-land model had been used. 
Results of similar magnitude were obtained with the Nimbus 2 model. Figure C-2 compares two 
maps of the albedo of the Arctic poleward of 70° N during the period July 1 to 1 5, 1969, as obtained 
by use of the cloud-land and snow models. 

Conversely, over ocean (fig. C-3) the cloud-land model would underestimate the albedo by about 
3 to 4 percent if used instead of the ocean model. A somewhat larger underestimate would have been 
obtained with the Nimbus 2 model. Thus, it can be stated from this comparison that in evaluation of 
the Nimbus 2 data the albedo of the Earth-atmosphere system was underestimated by about 2 to 3 
percent over the almost cloudfree subtropical ocean areas. 

Figure C-4 shows a comparison between albedos obtained with the Nimbus 2 and the cloud-land 
models. The former causes 1 to 3 percent higher albedos over most areas where diffuse reflectance is 
higher than 10 percent. Thus, it should be expected that over such areas or the Sahara the results 
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from Nimbus 2 are about 2 to 3 percent higher. In a global average, Nimbus 3 data resulted in values 
between 28 and 29 percent, while the Nimbus 2 evaluation resulted in values around 31 percent 
(Raschke and Bandeen, 1970). 

A comparison between albedo values of measurements taken over Greenland and adjacent ice- 
pack areas derived from the Nimbus 2 and the snow models is shown in figure C-5. There the Nimbus 2 
model overestimated the albedo by more than 1 5 percent. 

Therefore it appears from this discussion that in evaluations of the Nimbus 2 data, the albedo 
over dark, cloudfree ocean areas might have been underestimated, while over all brighter areas, in 
particular over the Arctic, it might have been considerably overestimated. Thus, albedos of more than 
80 percent over Greenland might be too high. From the discussion of the spectral albedo of the 
Earth-atmosphere system over snow fields, the lower albedos found with the snow model appear to 
be more realistic than those from Nimbus 2 data. It should, however, be mentioned that the snow 
model may underestimate considerably the increase of the albedo of the Earth-atmosphere system 
over wide snow fields, because it has been obtained from ground observations (Kondratiev, 1965). 

It may be concluded from these considerations that the models used here are still insufficient, espe- 
cially over snow and ice. Only careful airplane measurements of the flux density of reflected solar 
radiation at different angles of illumination over specifically selected areas will provide the means to 
derive more accurate models. 
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ALBEDO WITH NIMBUS 3 SNOW-MODEL (percent) 


Figure C-1.— Albedo values obtained with cloud-land and snow models from same measurements 
over Greenland and Antarctica. 47° < f < 51°. (f = Sun's zenith angle.) 






Figure C-2.— Albedos {in percent) of the Arctic region obtained with the same Nimbus 3 measurements, July 1 to 15, 

1969. (a) Cloud-land model. 
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ALBEDO WITH NIMBUS 3 CLOUD-LAND MODEL (percent) 



2 4 6 8 10 12 14 16 18 

ALBEDO WITH NIMBUS 3 OCEAN MODEL (percent) 


Figure C-3.— Albedo values obtained with cloud-land and ocean models from the same measurements (Nimbus 3, 
orbit 861, day 168) over cloudfree ocean areas. r d < 10 percent; T b > 273 K, channel 2; 16°<f<42°). Numbers 
by each plotted point are r d values in percent. 
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Appendix D 


GLOBAL MAPS OF RESULTS FROM MEASUREMENTS DURING THE 
PERIODS MAY 1 TO 15 AND OCTOBER 3 TO 17, 1969 


In this appendix are shown maps (figs. D-l to D-8) of the albedo of the Earth-atmosphere 
system, of the absorbed solar radiation, of the outgoing longwave radiation, and of the radiation 
balance of the Earth-atmosphere system calculated from Nimbus 3 measurements obtained during 
one semimonthly period in the spring and fall seasons (both hemispheres, 1969). Solid isolines 
in all maps refer to daily averages of the indicated quantities. 
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Figure D-1. -Albedo (percent) of the Earth-atmosphere system during the period May 1 to 15, 1969. Values are daily averages. 
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Figure D-3.-Solar radiation (calories per square centimeter per minute) absorbed in the Earth-atmosphere system during the period May 1 to 15, 1969. 

Values are daily averages. 
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Figure D-4.— Solar radiation (calories per square centimeter per minute) absorbed in the Earth-atmosphere system during the period October 3 to 

1969. Values are daily averages. 
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period October 3 to 17, 1969. Values are daily averages. 
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Figure D-7.— Radiation balance (calories per square centimeter per minute) of the Earth-atmosphere system during the period May 1 to 15, 1969. Values 

are daily averages. 
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Values are daily averages. 




Appendix E 


RADIANT EXITANCE W OF A BLACKBODY 
COMPUTED FROM THE STEFAN-BOLT2MANN LAW 


m 

W 

fl 

W 

f 

W 

■a 

W 

190 


215 

0.174 

240 

0.270 

265 

0.401 

191 

m 

216 

All 

241 

.274 

266 


192 

.110 

217 

.180 

242 

.279 

267 

.413 

193 

.113 

218 

.184 

243 

.283 

268 

.419 

194 

.115 

219 

.187 

244 

.289 

269 

.426 

195 

.118 

220 

.190 

245 

.293 

270 

.432 

196 

.120 

221 

.194 

246 

.298 

271 

.439 

197 

.122 

222 

.197 

247 

.303 

272 

.445 

198 

.125 

223 

.201 

248 

.308 

273 

.452 

199 

.128 

224 

.205 

249 

.313 

274 

.458 

200 

.130 

225 

.208 

250 

.318 

275 

.465 

201 

.133 

226 

.212 

251 

.323 

276 

All 

202 

.135 

227 

.216 

252 

.328 

277 

.479 

203 

.138 

228 

.220 

253 

.333 

278 

.486 

204 

.141 

229 

.224 

254 

.338 

279 

.493 

205 

.144 

230 

.228 

255 

.344 

280 

.500 

206 

.146 

231 

.232 

256 

.349 

281 

.507 

207 

.149 

232 

.236 

257 

.355 

282 

.514 

208 

.152 

233 

.240 

258 

.360 

283 

.522 

209 

.155 

234 

.244 

259 

.366 

284 

.529 

210 

.158 

235 

.248 

260 

.372 

285 

.536 

211 

.161 

236 

.252 

261 

.377 

286 

.544 

212 

.164 

237 

.257 

262 

.383 

287 

.552 

213 

.167 

238 

.261 

263 

.389 

288 

.559 

214 

.171 

239 

.265 

264 

.395 

289 

.567 


W = aT 4 cal cm -2 min 1 


where o = 8.130 544 94 X 10“ n cal cm -2 min -1 KT 4 . 
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